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Using Monte Carlo simulation method in the canonical ensemble, we have studied the commensurate- 
incommensurate transition in two-dimensional finite mixed clusters of Ar and Kr adsorbed on graphite basal 
plane at low temperatures. It has been demonstrated that the transition occurs when the argon concentration 
exceeds the value needed to cover the peripheries of the cluster. The incommensurate phase exhibits a simi¬ 
lar domain-wall structure as observed in pure krypton films at the densities exceeding the density of a perfect 
(\/3 X \/3)R30° commensurate phase, but the size of commensurate domains does not change much with 
the cluster size. When the argon concentration increases, the composition of domain walls changes while the 
commensurate domains are made of pure krypton. We have constructed a simple one-dimensional Frenkel- 
Kontorova-like model that yields the results being in a good qualitative agreement with the Monte Carlo results 
obtained for two-dimensional systems. 
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1. Introduction 


The commensurate-incommensurate (C-IC) transitions in adsorbed films have been experimentally 
observed in a variety of systems GHil and studied by theoretical methods and computer simu¬ 

lations i2^l3^ . 

In a large number of systems, the incommensurate phase forms a collection of finite commensurate 
regions (domains) separated by domain walls IHIl . The domain wall networks appear when the C phase 
is degenerated, i.e., when the atoms (or molecules) in the C phase occupy one of several equivalent lat¬ 
tices. The stability of IC phases with domain walls networks is considerably enhanced by large entropic 
effects The most mominent example of such a system is the krypton monolayer adsorbed on the 

graphite basal plane ISaldlll . At low temperatures, the adsorption of krypton leads to the formation of the 
three-fold degenerate commensurate (\/3 x \/3)R30° (\/3 x \/3)R30° structure, in which every krypton 
atom occupies one of every three hexagons of the graphite surface, has the lattice constant \/3 times the 
graphite lattice constant (ngr = 2.46 A) and is rotated with respect to the graphite lattice by 30 degrees. 
Upon the increase of pressure, the commensurate phase undergoes a transition into the incommensurate 
phase of higher density, in which the commensurate domains are sep arated by walls, and the atoms in 
neighboring domains are shifted to a different sublattice j4^|42-^|. In this case, the walls are called 
heavy and carry all the excess of density. At the densities just above the C-IC transition the domains are 
very large and their size gradually decreases when the density increases towards the monolayer comple¬ 
tion. This has been clearly demonstrated by a large scale computer simulation study (^. 

The C-IC transition has also been experimentally observed in the mixed Kr-Xe and Ar-Xe films on 
graphite li^liil . On the other hand, the structure of mixed Ar-Kr films on graphite has not been exper¬ 
imentally studied yet. The only available experimental results were published by Singleton and Halsey 
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nearly 60 years ago and demonstrated that mixed liquid-like films exhibit complete mixing, indepen¬ 
dently of the composition. The only information about the structure and properties of solid-like mixed 
Ar-Kr films on graphite have been obtained using computer simulation methods In particular, it 

has been demonstrated that already at submonolayer coverages and at sufficiently low temperatures, the 
addition of argon into krypton film triggers the C-IC transition, and leads to the formation of finite clus¬ 
ters in which small commensurate domains consisting of pure krypton are separated by heavy walls built 
of argon and kr 5 rpton js^l- The C-IC transition in one-component finite clusters of Lennard-Jones atoms 
adsorbed on graphite has also been studied by Houlrik et al. (i^l. They have used Monte Carlo method to 
study the effects of the corrugation potential on the structure of finite systems and have shown that for a 
given amplitude of the corrugation poteittial the presence of free surfaces enhanced the stabihty of the C 
structure. 

In this paper we study the mechanism of the C-IC transition of submonolayer mixed films made of 
argon and krypton using Monte Carlo simulation methods and appropriately modified one-dimensional 
(ID) Frenkel-Kontorova (FK) model The original FK model handles an infinite chain of atoms 

interacting via harmonic potential and subjected to a periodic (sinusoidal) external field at zero temper¬ 
ature. The FK model can be used to describe the basis features of C-IC transitions & At this point, we 
should mention that the FK model has been extended to two-dimensional systems Issll to mixtures 
( 5 ^, systems with disorder and has also been used to study finite chains 

In one of our recent works (iill . we studied the impurity driven commensurate-incommensurate tran¬ 
sitions in one-dimensional finite systems using the FK model. We assumed that the pure chain is commen¬ 
surate and discussed the effects of impurities located either only at one end or at two ends of the chain. 
It is shown that in both situations the C-IC transition occurs when the amplitude of the external field 
experienced by the impurity atoms falls into the region between the lower and upper threshold values. 
These limiting values of the amplitude depend upon the parameters characterizing the interaction be¬ 
tween the atoms in the main chain, the amplitude of external field acting on the main chain atoms and 
the interaction between the main chain and the impurities. The number of solitons (domain walls) in the 
IC structure were found to depend upon the chain length, and the parameters describing the interactions 
in the system. Our findings have been confirmed by Monte Carlo simulation in two-dimensional systems. 

Here, we also considered finite chains of Kr atoms with impurities (Ar atoms) located at the chain 
ends as well as in the chain interior. 

The paper is organized as follows. In the next section we presented the two-dimensional model and 
described the Monte Carlo method used. The results of Monte Carlo simulation are presented in section 3. 
Then, in section 4 we presented the modified ID FK model and demonstrated that it yields the results 
qualitatively very similar to those emerging from the two-dimensional Monte Carlo simulation. 


2. The two-dimensional model and simulation method 


We have considered strictly two-dimensional mixed submonolayer films consisting of Ar and Kr 
atoms. The interaction between the atoms is modelled via the truncated (12,6) Lennard-Jones potential 


Uij{r) = 



^ r max > 
^ > f^max > 


( 2 . 1 ) 


where r is the distance between a pair of atoms, the subscripts i and j mark the species Ar and Kr, and 
the potential is cut at the distance rmax = 3.0aij. The potential parameters for a pair of unlike atoms are 
given by the Lorentz-Berthelot relations: 

1 ,_ 

O'ylrJCr - - (O'Ar + O'Kr) and CArKr - (2.2) 

The external field due to the graphite substrate is assumed to be given by fi^lilll 


v[x, y) = - 14 ,/ {cos(qir) -i- cos(q 2 r) -i- cos[(qi - q 2 )r]}, 


(2.3) 
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Table 1. Lennard-Jones parameters for Ar and Kr used in this work 


ij 


ei,t,K 

Ar,Ar 

3.4 

120.0 

Kr,Kr 

3.6 

171.0 

Ar,Kr 

3.5 

142.83 


where qj and q 2 are the reciprocal lattice vectors of the graphite basal plane and the amplitudes 14 ,; (i = 
Ar,Kr) determine the potential barriers between adjacent minima. We assumed here that - 8.4 K 
and = 14.4 K jsil]. 

The values of parameters entering the potential <2.1) are given in table[l] The graphite lattice constant 
flgr = 2.46 A is assumed to be a unit of length and Eat is taken as a unit of energy. 

The Monte Carlo simulation is carried out in the canonical ensemble for systems of different 
(submonolayer) density, of different total number of atoms, N, and of different mole fraction of argon, 
xat, between 0 and O.SS. We used a standard Metropolis sampling and two types of moves are taken into 
account, i.e., the translation of a randomly chosen atom by a randomly chosen vector within the circle of 
radius dmax and the identity exchange. The magnitude of rfmax is updated every 1000 Monte Carlo steps 
(each Monte Carlo step consists of N attempts to move a single atom as well as N attempts to change 
the identity of randomly chosen atoms) in order to keep the acceptance rate equal to about 0.3. The 
equilibration and production runs consist of 5 ■ 10® and 5 ■ 10^ Monte Carlo steps, respectively. Standard 
periodic boundary conditions were used therein. 


3. The results of Monte Carlo study of Ar-Kr submonolayer films 

Since our study was performed for finite clusters at low temperatures it seems reasonable to begin 
with the estimation of the effects of free boundaries on the ground state properties of finite commensu¬ 
rate patches made of pure krjrpton. 

We performed Monte Carlo calculations at low temperatures using different starting configurations. 
At first, the rectangular simulation cells of the size Lx x Ly - L x L \/3l2 with L - 30, 60 and 120, were 
used. The number of Kr atoms in a perfect and fully covered {i/3 x \/3)R30° commensurate structure 
was equal to A/rt = L^I3. Then, we modified the simulation cell by making Lx larger (Lx - L + A), while 
keeping the number of Kr atoms unchanged, i.e., equal to A/rt = i^/3. In this way, two free boundaries 
running along the y-axis were created. Of course, the distance between these two boundaries, given by A, 
should be larger than the cutoff distance of the Lennard-Jones potential [equation <2.11 1 . The second series 
of calculations were done for finite hexagonal clusters of the commensurate phase built of a different 
number of krypton atoms. The simulation was carried out at reduced temperatures between 0.01 and 
0.1 and the energies obtained were extrapolated to T = 0. In this way, we estimated the ground state 
energies. 

In the case of a fuUy occupied rectangular box, we did not observe any finite size effects, and the 
ground state energy (per particle) of the commensurate phase was estimated to be equal to eo,c = -4.617 
(see figure [1). Similar calculations carried out for the systems with free boundaries along the y-axis led 
to the ground state energy equal to eo,f = -4.562 (see figure [1) for three different systems of the size 
(60 X 15 V3) (AIrt = 600), (60 x 45 n/S) (Art = 1800) and (90 x 301/3) (Art = 1200). 

In the ground state, the excess interfacial energy per unit length due to the presence of free bound¬ 
aries of the total length Lint can be calculated as follows: 

f^int — (fS’o.f — ! Aint i (3.1) 

where Lo.f = AiRpeoy and Lq.c = A/Rreo,c are the total energies of the systems with and without free bound¬ 
aries. For the three systems mentioned above, we obtained the same value of the excess interfacial energy 
at r = 0 equal to about 0.635. This indicates that finite size effects are negligibly small at very low tem¬ 
peratures. 
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Figure 1. The temperature changes of the energy (per particle) of fuUy covered (open circles) and with 
free boundaries (open diamonds) commensurate phase obtained using a rectangular slabs containing 
600,1200 and 1800 atoms and of different total length of the free Interface, = Ly Vs, with Ly - 30, 60 
and 90. The fiUed circles are the ground state energies of finite hexagonal patches plotted against 


On the other hand, the ground state energies of finite hexagonal clusters are bound to show large 
finite size effects. In such clusters, the fraction of atoms at the patch boundaries is proportional to 1/ \/lV. 
In order to extrapolate the results to the limit Afio- —" oo, we made the plot (given in figure[T) of the ground 
state energies (per atom) for the clusters of different size ehex,o(Mcr) against 1/ s/N and extrapolated the 
results to 1/ \/N - 0. In this limit, the boundary effects should vanish and the energy for ATkt ^ oo should 
be equal to eo.c- The obtained value is equal to -4.169 and agrees very well with the value of eo,c obtained 
for the fully covered surface. The above presented results demonstrate that our simulation method gives 
reliable results. 

It was demonstrated in our earlier work (illl . that finite submonolayer clusters made of the mixture of 
Ar and Kr atoms exhibit domain-wall structures, in which the commensurate domains made of kr 5 q)ton 
are separated by the domain walls of varying composition. Here, we considered such systems in a more 
systematic way, aiming at the determination of conditions necessary for the domain-waU structure to 
appear. 

The Monte Carlo simulations were carried out for a series of hexagonal clusters of a different total 
number of atoms equal to 343, 601, 931 and 1333 and of also different argon mole fraction between 0.08 
and 0.55. The calculations were performed over a rather wide range of reduced temperatures between 
0.01 and 0.45. Figure[2]shows the temperature changes of the potential energy (per particle) in the systems 
of different cluster size and of different composition. Only in the case of xat = 0.1 and the clusters with 
N up to 931, the energy smoothly changes with temperature, and the inspection of snapshots showed 
that aU krypton atoms form a commensurate patch while all argon atoms are located along the patch 
boundary [see figure[^(a)]. 

In order to dis ting uish the commensurate and incommensurate atoms, we used the following order 
parameter jilili^lsifl 

- cosiq^r) -i- cos(q 2 r) -i- cos[(q^ - q 2 )r], (3.2) 

and assumed that the atom is commensurate (incommensurate) when (p>0i(p^0). 

In aU other cases, the domain-wall structures were observed at sufficiently low temperatures. Fig¬ 
ure [^(b) shows the example of a snapshot recorded at T - 0.02 and x^r = 0.1 for the larger patch made 
of 1333 atoms. It can be clearly seen that unlike in the smaller cluster of Af = 931 atoms, the cluster with 
N = 1333 exhibits the domain-waU structure, and the walls are made nearly entirely of krypton atoms. 
The argon contributes to the walls only in the areas close to the patch boundaries. The patch boundaries 
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Figure 2. (Color online) The temperature changes of the potential energy (per particle) obtained for finite 
hexagonal patches of different size and of different argon mole fraction. The vertical blue lines mark the 
locations of transition between the mixed commensurate structure (stable at high temperatures) and the 
domain-wall structure (stable at low temperatures) in the largest clusters of N = 1333 atoms, obtained 
from annealing and freezing runs. 


are covered by a single layer of argon atoms. The structure of the 1C phase is qualitatively similar to the 
1C phase in pure krypton monolayers of the density exceeding that of a fully filled C structure. However, 
the mechanism leading to the commensurate-incommensurate transition is quite different in both situa¬ 
tions. In the case of a fully filled monolayer, the density excess over the value corresponding to a perfect 
C phase triggers the formation of domain walls. Upon the increase of density, the size of commensurate 
domains gradually decreases (^|^. Ultimately, a dense 1C phase without commensurate domains is 
formed close to the monolayer completion. 

In finite mixed clusters the atoms close to the boundaries are strained and the magnitude of this strain 
depends on the number of argon atoms in the system. It is shown that the excess interfacial energy in 
pure krypton clusters is positive and it induces a force normal to the interface and pointing towards the 
patch interior. Evidently, this force is too weak to lead to the C-IC transition in the pure krypton clusters. 
The presence of a layer of argon atoms at the boundaries causes the interfacial energy to be considerably 
higher and hence induces a considerably larger force, which is sufficient to compress the patch and to 
trigger the C-IC transition. The number of argon atoms necessary to cover the entire patch boundary is 
proportional to the total length of the interface, which in turn is proportional to For a small and 
fixed argon mole fraction in the system, the concentration of argon along the interface increases with the 
cluster size. Figure[^(a) shows that in the case of TV = 931 and xat = 0.1, argon atoms hardly cover the 
entire interface. On the other hand, when N - 1333 [figure[3](b)], the entire interface is densely covered 
by argon and there is some excess of argon atoms which try to penetrate the patch and trigger the de¬ 
velopment of domain walls. When the argon concentration increases, the walls become mixed and for 
a large argon mole fraction are predominantly built of argon (see figure |4). We studied the changes of 
the waU composition resulting from the changes of argon concentration in the system for the clusters of 
different size with N - 601, 931 and 1333. We did not take into account the incommensurate atoms with 
less than 5 nearest neighbors. Such atoms are located along the patch boundaries, and the inspection of 
snapshots (see figure[3]and|4) allowed us to assume that there is a single atomic layer of incommensurate 
atoms at the patch boundaries. Using such a procedure we were able to monitor the changes of the num¬ 
ber of argon and kr 5 rpton atoms incorporated within the walls and to estimate the argon concentration 
within the waU and along the patch boundaries. We performed a series of Monte Carlo simulations over a 
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Figure 3. (Color online) The examples of configurations for the finite patches with xpj- = 0.1 at T = 0.02 
consisting of 931 (part a) and 1333 (part b) atoms. Blue and red circles mark krypton and argon, respec¬ 
tively. The open circles correspond to the commensurate atoms with the order parameter (piri) & 0, while 
the filled circles mark the incommensurate atoms with (piri) < 0. 

certain temperature range between 0.01 and 0.3 for different cluster sizes and different argon mole frac¬ 
tion in the film. The examples of temperature changes of the wall composition are given in figure [S] (a), 
which shows the numbers of argon and krypton atoms forming the walls in the clusters of 1333 atoms 
and different argon mole fraction. The results demonstrate that for xat = 0.1, the patch is commensu¬ 
rate at sufficiently high temperatures, and undergoes a transition to the incommensurate structure at 
T = 0.14 + 0.02. The transition seems to be of the first order as suggested by rather large hysteresis along 
the freezing and annealing runs. In the case of xat = 0.2, the transition occurs at higher temperature 
of about T - 0.20 + 0.02 and is rounded. For the patches with a higher argon concentration, the incom¬ 
mensurate structure occurs over the entire range of temperatures studied. It is also evident that upon 
the increase of argon concentration in the film, the composition of walls changes from krypton rich to 
argon rich. Having the temperature changes of the wall composition we could estimate the ground state 
behavior of the systems studied. 

The results are given in figure (b), which shows the changes of the argon mole fraction within 
the walls and at the patch boundaries versus the total argon mole fraction in the system. These results 
demonstrate that when the cluster size increases, the onset of the C-IC transition occurs at a lower argon 
concentration in the film. In the cluster with N = 601, it occurs when x^r ~ 0.125, In a larger cluster 
consisting of 931 atoms it starts when xat is slightly higher than 0.1, while In the largest cluster of TV = 



Figure 4. (Color online) The example of configuration for the cluster with N = 1333 atoms and xat = 0.3 
at r = 0.06. The labeling of atoms is the same as in figure[^ 
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Figure 5. (Color online) (part a) The temperature changes of the numbers of argon and krypton atoms 
within the walls in finite clusters of N = 1333 for the systems with different argon mole fraction equal 
0.1 (circles), 0.2 (squares) 0.3 (diamonds) and 0.4 (triangles). Black and red symbols mark the numbers 
of argon and krypton atoms, respectively. (Part b). The mole fraction of argon within the wall (black 
symbols) and the mole fraction of argon atoms at the patch boundaries (red symbols) against the total 
mole fraction of argon atoms obtained for clusters of different total number of atoms (given in the figure). 
The lines serve only as a guidance. 


1333 the transition is found already when Xat ~ 0.09. Of course, when the argon mole fraction is lower 
than the limiting value necessary to trigger the C-IC transition, all argon atoms are located along the patch 
boundaries. The argon mole fraction within the walls increases linearly with the argon mole fraction in 
the film and for sufficiently large argon concentration the walls are entirely formed by argon. A further 
increase of the argon mole fraction leads to a gradual decrease of the size of kr 5 rpton commensurate 
domains j^. 

The inspection of configurations showed that when the walls are predominantly made of krypton 
atoms (at low argon concentrations) they are wider than the walls formed in the films with high argon 
concentration. The domain walls made of argon atoms only were observed to usually consist of just two 
rows of argon atoms and can be classified as heavy walls. The neighboring domains are shifted by the 
displacement vector of the length equal to Ugr (ill]. The results obtained for different cluster sizes showed 
that the domains are rather small, and usually contain up to 10 commensurate krypton atoms along the 
line joining opposite walls. This is quite different from the already mentioned pure kr 5 q)ton 1C phase, in 
which the domains are very large close to the C-IC transition and their size gradually decreases when the 
density becomes higher. 

We also studied semi-infinite systems consisting of the rectangular commensurate domain of the size 
Lx X Ly with periodic boundary conditions applied along the x and y-axes, but with free interfaces run¬ 
ning along the y-axis. The simulation cell was of the size Lx-G0 + A, with A = 30, and Ly - 30 \/3/2. 

Starting from a perfect C structure, we estimated the ground state energies for three different config¬ 
urations with the fixed numbers of Ar and Kr atoms equal to N^r - 60 and - 540. The simulations 
were carried out at low temperatures, between 0.01 and 0.1, using two types of algorithms, with and 
without the identity exchange attempts. 

Three different systems were considered. In the first (i), two rows of krypton atoms adjacent to each 
boundary running along the y-axis were replaced by argon atoms. The second system (ii) contained only 
single rows of argon atoms at both free boundaries and the rest of Ar atoms formed a single wall consist¬ 
ing of two rows of atoms inside the C patch at a certain distance ru, from one of the interfaces. The third 
system (iii) also consisted of single rows of argon atoms along the two free boundaries, but the rest of Ar 
atoms formed two walls, each made of one row of atoms located at the distance ru, from the closer free 
interface. The calculations showed that in the systems (ii) and (iii) energy becomes independent of ru, as 
soon as is larger than the assumed cutoff distance rmax- The results given below were obtained for 
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fw— 10 . 

The simulation without identity exchange attempts demonstrated that the system (ii) has a lower 
ground state energy (per particle) (ei -4.450) than the systems (i) {ea ~ -4.425) and (hi) (eui -4.398). 
Moreover, the entire krypton patch in the system (i) was observed to form the C structure at low tem¬ 
peratures, while the vast majority of argon atoms were displaced to incommensurate positions due to a 
certain contraction resulting from the presence of free interfaces. In the case of system (ii), argon atoms 
were observed to form an incommensurate wall with the displacement vector between the commensu¬ 
rate domains at both sides of the wall corresponding to a heavy wall. The structure of system (hi) was 
found to be commensurate, just the same as in the case of system (i) Another Monte Carlo run, in which 
the identity exchange attempted was taken into account, performed for the system (i) demonstrated that 
already at very low temperatures the system structure spontaneously changes to that corresponding to 
the system (ii). Moreover, the two single rows of argon atoms inside the patch of the system (hi) were 
observed to merge into a single heavy wall consisting of two rows of argon atoms. We also performed 
a run using the algorithm involving the identity exchange attempts and the staring conhguration with 
three rows of argon atoms at each free boundary. In this case, we observed the formation of two heavy 
waUs running along the y-axis and made of two rows of argon atoms again. 

We calculated the ground state energies of perfectly commensurate (iic.fc) and relaxed (iirei.fc) struc¬ 
tures for the systems k =(i), (ii) and (ih) as well as the energy gain (per the unit length of free interfaces) 
due to the relaxation 

®int,fc — (^^rel.fc “ ^C,fc) ^^int ■ (3.3) 

The energies (per atom) of perfectly ordered commensurate structures are very similar and equal to 
ec,i = -4.3748, ec ,2 = -4.3699 and ec ,3 = -4.36707. On the other hand, the energy gains due to relaxation 
are considerably different and equal to eint.i ~ -1.16, eint ,2 ~ -1.85 and eint ,3 ~ -0.76. 



Figure 6. (Color online) The example of configuration for the rectangular slab of the size 60 x 30 V3 with 
xat = 0.2 at r = 0.03. Blue and red circles mark krypton and argon respectively. Circles with thin lines 
correspond to commensurate atoms with the order parameter (^(r;) s 0, while circles with thick lines are 
for incommensurate atoms with 0(r/) < 0. 

The above presented results demonstrate that the formation of heavy walls consisting of two rows of 
argon atoms inside the krypton patch leads to the most stable structure. It should be also noted that sim¬ 
ulation for the systems with larger argon concentrations lead to the development of networks of heavy 
walls like that given in figure [6] rather than to a series of parallel walls running along the y-axis. This 
demonstrates that wall crossings contribute to the system stability and suggests that the waU crossing 
energy is negative. Of course, at finite temperatures the network of hexagonal walls is also stabilized by 
entropic effects, in particular, by the so-called breathing entropy 1^1^16611 . 

Concluding this section, we should mention that from our earlier study Isill it follows that upon the 
increase of temperature the domain walls structure disappears and the mixed commensurate phase is 
formed at the temperatures below the melting point. 
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4. A simple one-dimensional Frenkel-Kontorova model 

In this section, we propose a simple modified one-dimensional Frenkel-Kontorova model dHlil] that 
leads to quahtatively similar restilts as those discussed in the previous section. 

Taking into account that the domain-wall structures that appear for submonolayer coverages are 
stable only at very low temperatures, we consider finite chains at zero temperature. 

In general, the energy of the chain consisting of N atoms of two different species A and B can be 
written as 

^ fw-i 2 iV 1 

Y. [l-cos(2:?rj[:i/^z)] k (4.1) 

^ [ i=l i=l J 

In the above, we assumed that the interaction between a pair of nearest neighbors is harmonic and 
characterized by the force constant Kjj+i - K(A,A), K{A,B] or kr(B,B) and the equilibrium distance 
bij+i = h(A,A), h(A,B) or h(B,B), depending on the composition of the pair. The second sum in equa¬ 
tion <4.11 represents the contribution to the potential energy due to the periodic substrate field, with the 
distance between adjacent minima equal to a, and the amplitude of the external field Vi is equal to Va or 
Vb- 

In order to find the equilibrium configuration of such a mixed chain of the length N, one needs to 
specify the numbers of atoms A and B and their positions along the chain, and then to minimize the 
energy with respect to the set of {x,}. 

Here, we considered a few situations that are supposed to mimic the behavior of the previously dis¬ 
cussed two-dimensional systems. We assumed that the component A is Kr-like, while the component B is 
Ar-like. The equilibrium configuration of the chain made of atoms A only should correspond to the com¬ 
mensurate structure. On the other hand, the chain made of only B atoms should lead to an incommensu¬ 
rate floating structure. Assuming that the elastic constant K{A,A] is the unit of energy and the substrate 
lattice constant a is the unit of length, we found out that the model with the parameters Va - 0.0055 and 
b(A,A) - 1.9 describes reasonably well the commensurate structure. Note that by taking b(A,A) - 1.9 we 
assumed that in the commensurate structure every second potential weU is occupied by an atom and that 
the misfit itiaa - b{A,A] - 2 = -0.1 is negative. This value of misfit is larger than In the case of Kr/Graphlte 
system, for which it is equal to about -0.06 m, but assuming the value of the misfit equal to -0.06, 
the FK model leads to the stable commensurate phase even for unreasonably low amplitudes of the sur¬ 
face field. To model Ar atoms we assumed that Vb - 0.003, h(B,B) = 1.8 and kr(B,B) = 0.9 The parameters 
i<r(A,B) and h(A,B) are assumed to be given by the Lorentz-Berthelot mixing rules, so that 

is:(A,B) = \/li:(A,A)^:(B,B) and h(A,B) = 0.5 [h(A,A) -t h(B,B)]. (4.2) 

We considered the chains that begin and terminate with atoms B. This assumption arises directly 
from the observation that in two-dimensional systems, the domain-waU structures appearing only after 
the peripheries of the cluster made of Kr atoms are covered with a filled single layer of Ar atoms. The 
situation in which argon atoms are located only at the ends of the chain corresponds to the model labeled 
as Model I. We considered two versions of this model assuming that there is only one argon atom at each 
chain end (Model II) and that each end is decorated by two argon atoms (Model 12). Then, in the Model 
II, we assumed that at a certain position after the atom Ni, i.e., at the position Ali -i-1, a mixed wall of 
the thickness equal to k atoms appears. The version of the model with a given value of k is named Model 
IKfc.ATi). 

In order to consider mixed walls, we assumed that within the wall, the parameters representing the 
elastic constant, K^, the equilibrium distance, bw, and the amplitude of the substrate potential, Vw, are 
all dependent on the waU composition measured by the mole fraction of component B, xb, and are given 
by: 


= XBii:(B,B) + a- xb]K{A,A], 

(4.3) 

h^ = XBh(B,B) = (l-XB)h(A,A) 

(4.4) 

Vu, = XbVb+11-Xb)Va. 

(4.5) 
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The pairs (ATi, ATi + 1) and {Ni + k,Ni + k+ 1), i.e., at both sides of the wall, consist of one A atom and 
one wall atom, so that the respective elastic constants and equilibrium distances are given by 

r(A,m) = XB^:(A,B) + (l-XB)^:(A,A) and h(A, m) = XBh(A,B) = (1 - jCB)h(A,A). (4.6) 

The above assumption takes into account the fact that when the wall is made only of atoms B, the inter¬ 
action between the pairs {Ni,Ni + 1) and (Afi -i- fc, iVi -i- A:-i-1) is the same as between the pairs (A,B) and 
{N- 1, AT), which are equal to ki(A,B) and b(A,B), respectively. 

Of course, the above assumed additivity relations are too simple to properly describe the effects of 
waU composition on the behavior of adsorbed films, but we wUl demonstrate that even such a simple 
model leads to the results qualitatively very similar to those obtained for two-dimensional films. 

Then, we introduced the displacements w/ = Xj / b - 2 /, with b being the equilibrium distance for the 
pair i and i + 1, and rewrote equation <4.1) in the form 



(4.7) 


where the energy is expressed in units ofkr(A,A)a^, and.^i,/+i = Ki^i+il[K{A,A)a^] and 0/ = Vil[K{A,A]a^] 
and the misfits m, are defined as triij+i = b{a,p)la-2 with {a,p) - (A,A), (A,B), (A, w) or (w,w). 

In the case of Model 1, the energy given by equation (4.7) depends on the number of argon atoms at 
the ends of the chain, while in the Models II(fc, Ni) it depends on both the thickness of the wall (fc) and its 
position (Ail). 

One should note that the systems considered in the previous section were studied using Monte Carlo 
method allowing for the exchange of atoms identity, so that the distances between the patch boundary 
and the walls and the waU-wall distances inside the patch correspond to the equilibrium states. In the 
present one-dimensional model, the location of the wall is predetermined by iVi, and one needs to min¬ 
imize the energy with respect to all displacements {u,} as well as with respect to Ni in order to find the 
equilibrium state of the chain. On the other hand, the wall thickness and composition are controlled by 
the argon mole fraction in the chain. The first series of calculations is aimed at determining whether the 
assumed model predicts the formation of solitons when B atoms are located only at the chain bound¬ 
aries, i.e., when Xb = 0.0 within the wall. The calculations carried out for different chain lengths, with 
N ranging from 11 to 101, and using Model II always led to a single commensurate domain, i.e., to the 
same qualitative results as for pure A chain. The only effect of the B atoms located at the chain ends 
was an increase of atomic displacements near the chain ends (see figure [7). This result agrees very well 
with two-dimensional Monte Carlo simulations carried out for low xat, up to the value for which the 



50 


100 


X. 


Figure 7. (Color online) Atomic displacements versus atomic positions for finite chains consisting of 41 
atoms. Black circles and line correspond to pure chain of B atoms, whUe red and blue circles and lines 
correspond to the chains decorated by one and two atoms A at each end of the chain, respectively. 
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Figure 8. (Color online) The main figure shows atomic displacements versus atomic positions for finite 
chains consisting of 41 atoms with the ends decorated with one A atom and containing one A atom within 
the chain at the position Ni +l (given in the legend and marked by a filled circle). The inset shows the 
changes of the system energy with the position of atom A within the chain. 


entire boundary of commensurate Kr patch was covered by a single layer of Ar atoms. However, when 
the Model 12 was used, the introduced perturbation appeared to be sufficiently high to lead to the forma¬ 
tion of solitons, i.e., the domain-wall structure (see figure [7). From the results of the previous section it 
follows that already at very low temperatures, the Monte Carlo simulation, starting with the configura¬ 
tion containing two rows of argon atoms at the boundaries, leads to the formation of a wall consisting of 
two rows of Ar atoms inside the Kr patch, and the boundaries are decorated by only a single layer of Ar 
atoms. Therefore, it is quite likely that the structure with two B atoms at each end of the chain is not the 
equilibrium state for the chain with four B atoms. 

Then, we considered the systems described by Models ll(fc, Afi) with the waU of a thickness k located 
inside the chain after the atom Ni. At first, we assumed that the chain length is fixed (AT = 41) and studied 
the behavior of the model when the concentration of B atoms within the waU changes. 

This situation mimics finite submonolayer systems with the gradually increasing concentration of 
argon. In the two-dimensional systems, the number of atoms is proportional to N^. Our choice of N 
corresponds to clusters with the number of atoms equal to about 1500 + 100 i.e., not much larger than 
those with 1333 atom that were studied by Monte Carlo simulations. 

The calculations carried out within the Model 11(1, Aii) for the wall consisting entirely of argon (xb = 
1.0) demonstrated that when 7Vi ^ 6, there develops a single soliton located behind the wall (see figure[8). 
The energy of the system acquires a minimum value when Ni = 6 (see the inset a to figure[^. 

When Ail > 6, a single atom B inside the chain does not lead to the formation of a waU. Only a rather 
smaU perturbation appears close to the impurity atom (figure|^. Evidently, the commensurate stripe be¬ 
tween the boundary and the atom B is already wide enough to sustain the stress induced by the impurity. 

In the case of the Model 11(2, Ni) with N - 41, we found out that for Xb up to 0.04, the commensurate 
structure is a stable state for all possible values of Ni. However, already for Xb = 0.05, the system under¬ 
goes a transition to the domain-wall structure. Figure |9] shows the atomic displacements obtained for a 
series of systems with Xb between 0.05 and 1.0. Thus, the results indicate that even when the perturbation 
introduced by a smaU concentration of component B within the waU is weak, it leads to the development 
of solitons. Two types of behavior can be singled out. For xb e (0.05,0.6), there are three solitons and with 
the increase of xb, the B atoms gradually approach the position of the first soliton. In the second regime, 
only two sohtons appear and the B atoms are located right within the waU. 

Of course, the results obtained for a high concentration of atoms B in the wall, do not correspond 
weU to the 2D systems studied using Monte Carlo simulations. For a large concentration of atoms B, one 
would expect the formation of several mixed walls with lower values of xb rather than the appearance 
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Figure 9. (Color online) Atomic displacements versus atomic positions for finite chains consisting of 41 
atoms with the ends decorated with one A atom and containing two A atoms within the chain at the 
positions marked by black circles obtained for different compositions of the wall inside the chain and 
corresponding to the location of the wall that minimizes the system energy. 


of a single wall with a high value of xb- Nevertheless, the model presented catches the most important 
features of the C-IC transition in finite systems driven by the presence of impurities. 


5. Summary 

Using Monte Carlo simulation method we studied the C-IC transition which occurs at low tempera¬ 
tures in finite clusters of Kr atoms contaminated with argon. It is demonstrated that pure krypton clusters 
form only the commensurate structure. This is consistent with experimental data clearly demonstrating 
that the commensurate structure is stable as long as the density does not exceed the density of a fully cov¬ 
ered commensurate phase Idlll . Also, the clusters retain a commensurate structure when small amounts 
of argon are added to the system. It is established that argon atoms are predominantly located along the 
cluster boundaries and the C structure is stable up to the point when the entire cluster boundary is cov¬ 
ered with a single layer of argon atoms. As soon as the argon concentration becomes higher than this 
threshold value, the system undergoes the C-IC transition leading to the formation of the domain-wall 
structure. The walls separating the commensurate domains are found to be heavy and their composition 
to change with the argon concentration. In particular, we found out that the argon concentration within 
the waU changes linearly with the total argon mole fraction in the system. Thus, at a certain upper limit 
of the argon concentration in the film, the waUs are made of argon only, while the cluster boundaries 
are still covered with a single atomic layer of argon. A further increase of argon concentration causes 
the waUs made of argon and the argon layer along the patch boundaries to become thicker. It should be 
noted that the size of commensurate domains does not depend much upon the argon concentration and 
the domains are made of pure krypton. This demonstrates that the mixing does not occur in the C phase. 
On the other hand, the mixing does occur within the incommensurate walls over a rather wide range of 
argon concentration. 

The above described behavior can be attributed to the rather large changes in the interaction energies 
between different pairs (Kr-Kr, Ar-Kr and Ar-Ar) when the system undergoes the C-IC transition. In the 
C phase, the distance between neighboring atoms is determined by the structure of the graphite lattice. 
A considerably smaller misfit of kr 5 rpton compared to argon, causes the system to maximize the number 
of Kr-Kr pairs and to avoid the appearance of Ar-Kr and Ar-Ar pairs. In the incommensurate state, i.e., 
within the waUs, the energy of interaction between Ar-Kr and Ar-Ar pairs is considerably higher since 
the distances between nearest neighbors within the heavy waUs are considerably lower and closer to 
the locations of the LJ potential minima. This effect enhances the tendency toward mixing. Besides, the 
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stability of walls is further enhanced by rather large entropic effects. 

We also showed that similar C-IC transition emerges from a simple one-dimensional Frenkel- 
Kontorowa model applied to finite mixed chains. It is demonstrated that for appropriately chosen param¬ 
eters of the model, its behavior is qualitatively very similar to that obtained from Monte Carlo simulations 
of two-dimensional systems. 
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MexaHiSM (|>opMyBaHH}i crpyKTypi/i aomchhoY ctihki/i b 
cy6MOHOLijapoBi/ix n,niBKax Ar-Kr Ha rpa(|)iTi 

A. naipi/iKieB, C. CoKo.iioBCbKi/iM 

Biflflia MOfleaiOBaHHfl (()i3MK0-xiMi4Hnx npopeciB, xiMiHHi/iit cjtaKy/ibTer, yHiBepcnrer iMeni Mapii 
CK/iOflOBCbKOi-KKDpi, 20031 f\K)6n\H, noabipa 

BuKopMCTOByKDHi/i MBTOfl MOfleaKDBaHHti MoHTe-Kapao B KaHOHitHOMy aHcaM6ai, mu ppcn\py\nv\ nepexifl Bifl 
cniBBMMipHOi flo HecniBBMMipHOi (jjasM y flBOBi/iMipHMx CKiHHeHi/ix 3MiLiiaHnx KaacTBpax Ar i Kr, aflcop6oBaHnx 
Ha rpac|)iTHiCi Oaai/iCHiit naoipHHi npw Huabxnx reMneparypax. flOKaaaHO, ipo nepexifl BiflOyBaetbcti, xoai/i koh- 
peHTpapia aproHy nepeBHipye 3HaHeHHa, Heo6xiflHi flat! noKpnrra nepi/ic|)epm Kaacrepa. HecniBBi/iMipHa c|)a3a 
fleMOHCrpye nofli6Hy crpyKtypy flOMeHHOi ctIhkh flO roi, axa 6yaa cnocrepexteHa y naiBxax nncToro Kpnnro- 
Hy npia rycTnnax, ipo nepeBwipyKDTb ifleaabHy (\/3x \/3)i?30° cniBBWMipHy (fiaay, aae po3Mip cniBBnMipHwx 
flOMBHiB 3HaHHO He 3MiHioeTbca 3i 3MiHOK) po3Mipy KaacTBpa. Koan KOHpeHTpapia aproHy apocrae, cxaafl flo- 
MeHHHX CriHOK 3MiHK)eTbCa, B TOH HBC SK CtliBBHMipHi flOMCHH CKtiaflaiOTbCa 3 HHCTOtO KpHHTOHy. Mh no6yfly- 
Baaw npocry OflHOBi/iMipny MOfleab (PpeHKeaa-KoHTopoBOi, axa npHBOflWTb flO peayabtariB, ipo flo6pe axiCHO 
yaroflxyKDTbca 3 peayabraraMH Monre-Kapao, orpHManuMH flaa flBOBHMipni/ix ciacreM. 
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K/homobI cnoBa: aficop6i4ift cyMiine^i, nepexofti/i evj cniBBi/iMipnoiflo HecniBBi^MipHOicjiasi^, MOflentOBaHHB 
MoHTe Kapjio, CKiHneHi cmctbmi/i, MOflejib 0peHKeji^-KoHTopoBOi 




